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In the last decade, the use of an ensemble of radio pulsars to constrain the characteristic 
strain caused by a stochastic gravitational wave background has advanced the cause of 
detection of very low frequency gravitational waves significantly. This electromagnetic 
means of gravitational wave detection, called Pulsar Timing Array (PTA), is reviewed 
in this article. The principle of operation of PTA, the current operating PTA and their 
status is presented along-with a discussion of the main challenges in the detection of 
gravitational waves using PTA. 
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1. Introduction 

A distinguishing feature between the Newtonian gravity and General relativity (GR) 
is the existence of gravitational waves, which are implied by any relativistic theory 
of gravitation. Einstein proposed the notion of gravitational wave (GW) formally, 
soon after formulating GRP^ As these waves can propagate without any significant 
absorption, they carry information about large gravitating masses, often exclusively, 
and complementary to that provided by electromagnetic waves. Such information 
is important for the development of gravity theories and by extension, for a funda- 
mental understanding of the origin and evolution of our universe. 

GW astronomy has seen significant development in the last 4 decades since 
the first instruments were conceived!^ Sensitive terrestrial instruments, such as 
Laser Interferometer Gravitational- Wave Observatory (LIGO) and Virgo, are op- 
erational and space-based instruments, such as New Gravitational Wave Observa- 
tory (NGO/eLISA^and ASTRO-GW, 13 are proposed in the future to detect these 
waves. International Pulsar Timing Array (IPTA) is an ambitious collaborative in- 
ternational effort of the entire radio pulsar community to detect gravitational waves 
using pulsed electromagnetic radiations from pulsars. While LIGO (and advanced 
LIGO) is sensitive to a frequency range between 10 Hz to 100 KHz, NGO/eLISA 
would be sensitive between 0.1 Hz to few /zHz. IPTA and its constituent Pulsar 
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Timing Arrays (PTAs) , owing to their large arm length, probe GW in the nano-Hz 
range (very low frequency bancP"^- 300 pHz to 100 nHqf|.) and are complementary 
with the terrestrial and space based GW detectors. A review of the current status 
of GW detection with the PTAs is presented in this paper. 

Pulsars are highly magnetized rotating neutron stars emitting a train of narrow 
radio pulses. They are known for their remarkable stability of period, which repre- 
sents the rotation rate of the star. The large rotational kinetic energy reservoir (10 46 
ergs) of this 20-km diameter star, with typical mass 1.4 M , is responsible for its 
remarkable clock like behaviour. As clocks play an important role in characterizing 
gravity theories, radio pulsars have been used to test these theories to a remarkable 
precision.^ The first detection, albeit an indirect one, of the existence of GW was 
also provided by radio pulsars, soon after the discovery of the first double neutron 
star (DNS) system, PSR 61913+16.™ As of today, 9 DNS systems are knowrP 1 
apart from other binary systems involving radio pulsars, which continue to provide 
new tests for gravitational physics. PTA utilizes the clock-like mechanism of radio 
pulsars to construct a telescope for the detection of very low frequency GW. 

A brief historical account of PTAs is in order here. The idea to use pulsar tim- 
ing residuals to constrain the stochastic GW background (SGWB), considering root 
mean square (rms) timing residuals of about 150 /is over 5 years of data in long 
period pulsars such as PSRs B1919+21, B1933+16 and B2016+28, was first pro- 
posed by Detweiler^ in 1979 closely following a suggestion by Sazhirf^ in 1978. 
Hellings & Downs^ further analysed data from 4 slow pulsars using NASA's deep 
space network antennas to put a limit on SGWB to about 0.01 percent of the 
critical/closure density. A similar limit from B 1237+25 was arrived at by Romani 
& TaylorpS w ho also compared it with other methods discussed by Zimmerman 
& Hellings Millisecond pulsars (MSPs), pulsars with short period (ms) and low 
magnetic fields, provide much higher period stability and precision of measurements 
of the time-of-arrival (TOA) of the radio pulse. Hence, MSP timing started soon 
after the discovery of the first MSP B 1937+21^ and interesting limits from long 
term timing of MSPs were put by several groups the most stringent one com- 
ing from the latter authors. Work of these groups, particularly Kaspi et alP^ also 
highlighted challenges presented by an instrument, which uses pulsar clocks, to set 
a limit or detect GW, namely, errors due to ephemeris, Dispersion Measure (DM) 
variations and so on. Both Detweiler and Hellings & Dowd^ED had proposed 
using a pair or ensemble of pulsars as a means to improve sensitivity not only to 
SGWB, but also towards individual sources. This was further extended by Foster 
& Backer^ to propose a GW telescope using an array of pulsars. Currently, three 
such experiments implement a substantially developed form of this concept and are 
also sharing their data in a major international effort to address this issue. 



a http: / /astrod.wikispaces.com/filc/vicw/GW-classification.pdf 

b Dispersion Measure (peem -3 ) is defined as the integrated column density of electrons along the 
line of sight 
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Although such a celestial instrument was proposed more than 3 decades ago, ad- 
vances in instrumentation and our knowledge of systematic effects in pulsar timing, 
based on two decades of dedicated sustained timing of a large number of pulsars 
from Jodrell Bank and Parkes radio telescope and other smaller scale efforts, have 
accelerated efforts in the last ten years to realise such an instrument. A sustained 
international collaborative effort of the entire pulsar community, particularly in 
the last 5 years, appears to be bringing such an instrument excitingly close to an 
eventual detection. A review of these extensive efforts is presented in this paper. 

The plan of the paper is as follows. The conceptual and theoretical framework 
and the technique used in a PTA is described in Section [2] The challenges in the de- 
tection of the small amplitude of GW are discussed in Section [3] A PTA requires an 
ensemble of pulsars. The criteria for the source selection and the observing require- 
ments are presented in Section |U The current PTA initiatives are then presented in 
Section [5] followed by a review of the status of these efforts in Section [6] The con- 
cluding section presents some remarks on future desirable developments for PTAs. 

2. An Array of Pulsars as a Gravitational Wave Telescope 

A typical PTA attempts to detect correlated perturbation in flat space-time, caused 
by a passing GW, by measuring systematic deviations from an expected rotation 
model of a radio pulsar using the technique of pulsar timing. A brief discussion of 
these ideas is presented here. 

2.1. Pulsar Timing 

Pulsar timing is an observational technique used to determine accurately the rota- 
tion parameters of a radio pulsar. Here we briefly review the technique. The pulsed 
emission from a radio pulsar, averaged over several thousands of periods, is a unique 
signature of the pulsaJ^ and is called the average profile (AP) of the pulsar. The 
time-of-arrival (TO A) of a fiducial point on the AP is computed by cross-correlating 
the observed AP of a pulsar with a noise free template constructed from several ob- 
servations of the pulsar. The observed TOAs are then transformed to a proper time 
at the solar system barycentre (SSB) through a chain of transformations, involving 
pulsar position, proper motion, parallax, DM, relativistic clock and space correc- 
tions in the solar system, and corrected for arrival at the SSB. The corrected TOAs 
are then compared with those predicted from an assumed model of the pulsar clock 
to obtain timing residuals. The assumed model of the pulsar clock involves the 
rotation rate of the pulsar and its higher derivatives, and may involve, for binary 
pulsars, transformations similar to those outlined above from binary centre of mass 
to pulsar frame involving Keplerian, and post-Keplerian parameters, if the binary 
is relativistic. The timing residuals are then fitted in a least squares fit to improve 
the parameters of the assumed model. If the fitted model describes the pulsar clock 
well, one should obtain "white" or Gaussian distributed residuals at the end of 
the process representing only the measurement errors. Any systematics or larger 
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than expected errors indicate unmodeled parameters. The techniques and the stan- 
dard software, TEMPO(2), used for this analysis are described in Refs. [20H251 The 
perturbation produced in the metric by a passing gravitational wave is one such 
parameter. Thus, rms timing residual provides an upper limit on this perturbation 
and can be used to constrain the amplitude of GWPS 

2.2. Effect of GW on timing residuals 

A passing GW introduces small time varying changes in the proper separation 
between two points in space-time. Thus, the time of flight of photons of a received 
electromagnetic radiation varies introducing a shift in frequency, a fact first noted 
by Estabrook & WahlquistpS GW detectors, such as LIGO and LISA, use phase 
difference obtained with interferometric measurements, between the emitted and the 
reflected radiation from a coherent source, such as a laser, to measure these changes 
in the proper separation. On the other hand, PTA employs a natural coherent 
source, namely pulsars, by measuring its apparent pulse frequency. The fractional 
difference in the apparent pulse frequency, v{t), and the pulse frequency emitted 
by the pulsar, v$, is given by the following expression for a monochromatic GW, 
propagating in z direction in a coordinate system with earth at origirP^ 

( V ® - ^ = -(1 - cos(0)) coa(2<t>)[h(t) - h(t - die - dcos(6)/c)} (1) 
v 2 

where, h is the amplitude of the wave. The earth-pulsar sight-line is at an angle 
6 in the x-z plane and the principal polarisation of the wave makes an angle <f> 
with the x-axis. The left hand side of Eq. [1] is equivalent to the derivative of the 
timing residual, thus relating the timing residuals to the amplitude of the GW. In 
particular, the spectral density of the residuals is identical to that of the amplitude 
of the wave^ and can be obtained from the auto-correlation of the residuals. 

The duration used to obtain the expected mean square residuals in the above 
auto-correlation sets the frequency range of the GWs to which the PTA is sensitive. 
Typical durations quoted are in the range of 10 7 to 10 8 s, corresponding to arm 
lengths in excess of 10 15 km and frequencies of 1 to 100 nHz. This distinguishes 
PTA from other GW detectors, where the electromagnetic radiation travels a very 
small fraction of the GW wavelength. This range of frequency is also complementary 
to other GW detection experiments. 

2.3. Upper Limits on SGWB using PTA 

The SGWB can be specified in terms of three different quantities of interest, which 
are related to each other and can be estimated from the timing residual measure- 
ments. Much of the earlier work used the estimates of rms timing residuals from 
individual pulsars to obtain an upper limit on the energy density of the GW as a 
fraction of the critical density in an FRW (Friedmann- Robertson- Walker) universe 
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given b^SMi 

^(/) = |^^7 (2) 

where, Hq = 100 ho km s -1 is Hubble's constant with ho expressing the uncer- 
tainty in the Hubble constant (typically 0.5 < h < 0.65^). The dimensionless 
energy density, Q gw (f), has an uncertainty due to ho- The ensemble average, over 
all directions and sources, of Fourier component of the metric perturbation is de- 
fined as the spectral density, Sh(f)- It can be estimated from the spectral density 
of the timing residuals. Since S/j(f) has the dimensions of inverse frequency, a di- 
mensionless quantity, h c (i), characteristic of strain amplitude, is often used and is 
related to 5/j(f) and il gw (i) as below 

h*(f) = 2fS h (f) 

3H^n gw (f) (3) 

27T 2 / 2 

The characteristic strain is related to the spectra of the timing residuals^! 

h 2 c (f) = 127r 2 / 3 i?(/) (4) 

where, R(f), is the power spectrum of the timing residuals. The variance of arrival 
time residuals, er 2 , is related to R(f) by the following expression 

/>oo 

* 2 = / R(f)df (5) 



In general, there are many contributions, other than that due to GW, to the 
noise seen in the timing residuals (See Section [3]). Hence, the estimate arrived with 
Eqs. 0] and [5] can only provide an upper limit to the characteristic strain, h c . Thus, 
individual pulsars with small timing residuals and "white" residuals provide the best 
upper limits. Much of the earlier work on the detection of SGWB used this method 
to estimate an upper limit by averaging the measurements of cr 2 for a sample of 
pulsars, assuming each to be an independent measurement. 

The change in apparent frequency, given in Eq.[I] incorporates information about 
GW amplitudes at two distinct space and time coordinates. The first term incor- 
porates information about h at the time and place of reception (at Earth) and the 
second at the pulsar. The first term will be correlated for many pulsars, whereas 
the factors contributing to the timing noise (See [3]) of two pulsars are generally 
uncorrelated. Monitoring of timing residuals in this coherent fashion from a large 
sample of pulsars is what constitutes a PTA. The aim of PTAs is to detect this 
correlation as a function of angular separation between pairs of pulsars in the PTA. 

The form of this correlation, ((9 ab ) @, is given bjffSMEU 

(1 - COS(M) (1 ~ COsjdab)) , 1 (1- COS(M) 1 (Q) 

^ ab) 2 2 ' 6 2 3 



c The scaling factor used in this definition for f varies in the literature. We have adopted the 
definition from Ref. 1131 
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Fig. 1. Expected correlation between timing residuals for a pair of pulsars as a function of sepa- 
ration angled 

which is plotted in Fig. [1] The angular separation between the pair of pulsars 
is 9 a b- Timing observations of a sample of pulsars can be analysed to look for a 
correlation, which provides a constraint on the amplitude of the GW spectrum. 
The characteristic strain spectrum for different models of SGWB can be generically 
written aJ^l 

hc(f) = Af a (7) 

Thus, a limit can be obtained in principle by a fit to correlations in the timing 
residuals of several pairs of pulsars with varying angular separations.^ 

2.4. Models of SGWB and characteristic strain spectrum 

While ground and space-based detectors are sensitive to GW generated by sources, 
such as supernovae, spinning neutron stars, in-spirals of extreme mass binary sys- 
tems, coalescing binaries and binaries with neutron stars and black holes and super- 
massive black hole binaries (SMBH), the most probable sources for PTA are SMBH. 
Galaxy mergers during the history of universe can lead to the central SMBH in these 
galaxies to become gravitationally bound as a SMBH binary and generate GW dur- 
ing its late time in-spiral. A PTA can potentially detect GW from such SMBH with 
orbital periods of 10 8 s. 

A large number of SMBH are predicted by the theories of formation and evo- 
lution of SMBH. The most likely source for PTAs is the SGWB resulting from an 
incoherent addition of GW generated by a number of such SMBH systems! 29 ^ 34 ! The 
characteristics strain spectrum for such a background is given by Eq. [7l with the 
power law index, a = -2/3. 29 

Alternatively, PTAs are sensitive enough to detect SGWB due to other models. 
One such alternative is the quantum-mechanical generation of cosmological pertur- 
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bationsP^The power law index, a, in the Eq.[7]is -1 corresponding to the parameter 
ft = -2 characterising the power law evolution of the cosmological scale factor, a, 
with the cosmic time (See Eq. 2 of Ref . l35| . Another alternative is SGWB formed 
by GW emitted by a network of cosmic strings The index in this case is -7/6. 
Therefore, PTA results can put meaningful constraints on these models. 

3. Factors affecting PTA precision 

Although more than 2000 pulsars have been discovered till date, it is difficult to 
obtain an ideal ensemble of pulsars, which satisfy all the requirements as listed in 
Section [4721 to form a high precision PTA. Apart from the pulsars being typically 
weak radio sources, the pulsed emission from pulsar exhibits a variety of instabili- 
ties, which makes extracting extremely tiny systematic deviations, caused by GW 
perturbation of metric, very challenging. These factors are discussed in this section. 

3.1. Rotational instabilities 

Although pulsars broadly exhibit clock-like stability, their rotation is affected by 
small rotational irregularities. Many young pulsars (characteristic age, r c < 10 6 
yr show sudden increase in their rotation rates, known as glitches P^*^ The fre- 
quency change at the glitch and its post-glitch recovery complicates timing analy- 
sis. Moreover, these events occur randomly. Hence, such pulsars are excluded from 
PTAs. 

None of older pulsars (r c > 10 8 ) have been reported to glitch and these tend to 
have periods larger than few 100 ms. While such pulsars were previously used to put 
constraints on SGWBP^the longer observations required for these pulsars to form 
a stable AP and changes in AP due to a variety of effects discussed later preclude 
their inclusion in the sample of pulsars for the PTA. Another class of old pulsars 
not showing glitches are the millisecond pulsars (MSPs) , which are pulsars with low 
surface magnetic fields (B < 10 10 G). The short periods help in obtaining a stable 
AP in short observations. Hence, MSPs are ideal candidates for a PTA sample. 

Many MSPs show slow wander of rotation rate over intervals ranging from 
months to years. This phenomena, called timing noise, is also not predictable.^ 
Timing noise is "red" in character and is similar to the expected GW signature 
and presents a significant constraint Recently, quasi-periodic structures in tim- 
ing residuals have been reported suggesting an origin intrinsic to the neutron star^ 
(See Section l3~4|) . Hence, a systematic characterization of a much larger sample of 
MSPs, than being used in the present day PTAs, is required to choose the most 
ideal rotators. 



d estimated from measurements of periods and period derivatives assuming dipole radiation from 
a pulsar 
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Fig. 2. The TOA delay due to an error of 0.08 pccm" 3 in assumed DM. Each panel shows the 
pulse obtained at 613 MHz and 237 MHz respectively with the GMRT. The upper panel shows 
the averaged pulse at the two frequencies with an assumed DM of 35.665 pccm -3 , while the lower 
panel shows that with the estimated DM (35.740 pccm -3 ). Figure reproduced from Ref. 1441 

3.2. Signal-to-noise considerations 

The TOA precision depends on the radiometer noise (and signal-to-noise ratio; 
S/N)P^ Hence, it is important to include strong pulsars in a PTA. On the other 
hand, radio pulsars are weak radio sources (typical flux density 10 mJy; 1 mJy = 
10 -29 W/m 2 /Hz) and this restricts the number of usable pulsars. Nearby MSPs are 
usually chosen for most PTAs. While the higher sky background leads to higher 
radiometer noise at low radio frequencies, pulsars are stronger at these frequencies 
due to their steep spectra. Thus, the best trade-off is achieved at frequencies > 1 
GHz. 
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3.3. Propagation effects 

The deviation in pulsar clock are measured using the pulsed radio signal from the 
pulsar, which propagates through the inter-stellar medium (ISM). Radio emission 
is dispersed and scattered by the electrons in the ISM and encounters multi-path 
propagation. This smears the pulse, reduces the effective time resolution and S/N 
and also causes a frequency dependent delay in TOA. A combination of all these 
effects, due to the inhomogeneities in ISM, causes variability in pulse TOA, com- 
promising the precision of a PTA. An illustration of these effects is shown in Fig. [2] 
andEP 

As pulsars are fast moving objects,^ observations of pulsars encounter different 
cross sections of ISM on different epochs. Fluctuation of electron density over these 
sections causes the dispersive delay in TOA of pulsars to vary. DM variations of the 
order of 0.0001 pccm - have been established in many past 

studiesJ^KMin] which 
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correspond to variations in dispersive delay in TOA of the order of 212 ns at 1400 
MHz. In addition, the line of sight to a few PTA pulsars pass close to Sun at specific 
epochs in a year, causing an additional dispersive delay due to the solar wind of the 
order of 1 ^s within 7° of the SunPEH 

Appropriate corrections, estimated from observations at two or more frequen- 
cies within a few days of timing observations, for such variations are included in the 
timing model. The usual PTAs do not carry out such multi-frequency observations 
simultaneously, but the comparison of the results of such non-simultaneous obser- 
vations with those from simultaneous multi-frequency observations, in progress at 
the GMRT, shows that the estimation uncertainties are almost identical.^ 

The dispersion of the radio pulse also reduces the time resolution for the pul- 
sar observed with the traditional filterbank instruments due to a finite channel 
bandwidth. Hence, PTAs utilize instruments implementing coherent dedispersion 
in either hardware or software, which allow even better estimation of DM varia- 
tions. 

Multi-path propagation of the pulsar signal leads to variable broadening of AP 
(Fig. 13]), particularly at frequencies below 400 MHz. As this is not easily correctable, 
timing observations are usually carried out above 1 GHz, where this effect is neg- 
ligible for low DM pulsars due to steep (~ frequency dependence of scatter- 
broadening. 

Lastly, intensity variations, usually caused by Diffractive Inter-stellar scintilla- 
tion (DISS) and Refractive Inter-stellar scintillation (RISS), affect the TOA pre- 
cision. The former can be addressed by employing large bandwidths compared to 
scintle size, which is also desirable for higher S/N. 




Fig. 3. Scatter-broadcning at low radio frequencies. Both panels show coherently dedispersed AP 
of MSP PSR B1937+21. The upper panel shows the AP at 610 MHz, whereas the lower panel 
shows the AP at 235 MHz, which is clearly scatter-broadened. Figure reproduced from Ref. 1441 
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3.4. Intrinsic changes in pulsed emission 

The pulsed emission from a pulsar is the probe of its clock mechanism provided the 
pulse (AP) is stable in its shape and arrival phase. Although this was broadly estab- 
lished in the early days of pulsar astronomy,^ observations over last four decades 
have shown that many pulsars switch between two or more stable APsPThis phe- 
nomenon is called profile mode-change. Recent studies indicate that profile mode- 
changes may be related to/or accompanied by changes in rotational parameters of 
a pulsar. Intermittent pulsars, such as PSRs B1931+24, J1841-0500, J1832+0029, 
J1839+1506, exhibit extended periods (several hundred days) with no pulsed emis- 
sion and different spin-down rates for their ON and OFF epochsP^HH] j± sam pi e 
of 17 pulsars with profile changes corresponding to changes in spin-down rate was 
reported recently suggesting the cause to be switching between different magneto- 
spheric statesP^Lyne et alP^l argued that the timing noise in MSPs can be explained 
by slow profile changes and claimed that it may be possible to correct for timing 
noise. In any case, the effect of the magnetospheric processes on the TOA precision 
poses a challenge, quite distinct from the processes related to the internal structure 
of neutron star, as discussed in Section 13.11 and motivates a better understanding 
of these processes to account for their contribution to timing residuals. 

3.5. Errors in timing technique 

In addition to the above effects, the timing residuals may have contributions from 
the clock errors in the atomic clocks used for obtaining the topocentric TOA and 
errors in the solar system ephemeris. See Ref. [17] for a detailed discussion on these 
errors. Both the ephemeris and the atomic standards are being steadily improved. 

4. Requirements for a PTA experiment 

The effects discussed in the last section influence the choice of pulsars and the 
observing instrumentation best suited for a high precision PTA. These requirements 
are summarised below. 

4.1. Observing requirements 

The appropriate instrumentation and observing strategy required for a typical PTA 
are 

• A telescope or combination of telescopes with large collecting area 

• Wide-band receivers, preferably with a bandwidth > 1—2 GHz 

• Observing frequency between 1 to 3 GHz for the best trade-off between the in- 
fluence of sky-background, steep pulsar spectra and the propagation effects 

• Simultaneous or near-simultaneous observations at 2 or more frequencies to cali- 
brate DM changes, with one frequency preferably below 400 MHz 

• Full polarimetric observations 
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• High time resolution receivers, employing coherent dedispersion, to eliminate the 
dispersion smearing of the pulse. 

The frequent monitoring of a large sample of pulsars using the most sensitive tele- 
scopes requires a large amount of observing time on these, usually over-subscribed 
telescopes. PTAs in operations currently try to share this burden, which was also a 
prime motivation for forming the International Pulsar Timing Array (IPTA). Use 
of different observing system introduces sometimes epoch dependent variable fixed 
pipeline delays and clock errors. Thus, an important requirement for instrumenta- 
tion is to periodically calibrate such delays. 

4.2. Sample selection 

The requirements for the ensemble of pulsars comprising an effective PTA are as 
follows 

• Pulsars with high signal-to-noise ratio (S /N) pulsed emission 

• Pulsars with stable AP 

• Pulsars with narrow pulses 

• Field MSP or MSPs in wide binaries are preferred over tight binaries to avoid 
relativistic effects in timing 

• Pulsars on line of sights without unusual ISM structures 

• Nearby pulsars are preferred over more distant ones 

• Pulsars with exceptionally high rotational stability 

• A reasonable distribution of pulsars mimicking the pairwise angular separation 
required for C(^ah) given by Eq. [5] 

The first six requirements have important implications for the precision of de- 
tection of a GW signature and are the consequences of the timing technique itself. 
The last two requirements are dictated by SGWB detection technique. 

The short periods of MSPs combined with sharp pulses allows averaging several 
thousands of their pulses. MSPs also have very stable rotation and a more uniform 
distribution across the sky as they form a relatively local population. Hence, all 
PTA samples use MSPs as their ensemble. 

Although more than 150 MSPs are known, the above criteria reduces the list of 
suitable pulsars to about 80. While this ensemble is useful, it is not ideal as most 
of these show systematics in timing residuals, such as timing noise. Currently, only 
PSR J1713+0747 shows white residuals over 5 years, with post-fit rms residuals of 
better than 200 nsP^OHl More such pulsars are required to be discovered to make 
an effective PTA. 

5. Pulsar Timing Arrays 

There are three currently active pulsar timing arrays - the Parkes Pulsar Timing Ar- 
ray (PPTA), the North American Nanohertz Observatory for Gravitational Waves 
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(NANOGrav) and the European Pulsar Timing Array (EPTA). A brief description 
of each of these PTAs highlighting the instrumentation in use is presented in this 
section. The details can be found in Refs. l33"I[50l[57l 

5.1. The Parkes Pulsar Timing Array (PPTA) 

PPTA, which started operating in 2004, was initiated by R. N. Manchester and is 
a collaboration between Australia Telescope National facility (ATNF), Swinburne 
University of Technology and several other Australian and international institutions. 
It uses 64-m Parkes radio telescope to time 20 MSPs every 2 to 3 weeks. The 
Parkes radio telescope can observe MSPs below a declination of + 25° and the 
entire southern Galactic plane, allowing it to observe a large sample of Galactic 
MSPs. Observations are carried out using 20-cm, 10-cm and 50-cm receivers, the 
latter two frequencies coming from a dual-band coaxial receiver. All these are wide- 
band receivers. Wide-band backends, such as Parkes Digital Filterbank systems 
(PDFB1/PDFB2 /PDFB3 / PDFB4) and baseband recording systems with coherent 
dedispersion, such as Caltech-Parkes- Swinburne-Recorder (CPSR2) and APSR are 
employed for the observations. See Ref. [5D] and references therein for details of 
the observing system, instrumentation used, data reduction software and off-line 
processing pipeline. 

5.2. The North American Nanohertz Observatory for 
Gravitational Waves (NANOGrav) 

NANOgrav is a collaboration between National Radio Astronomy Observatories 
(NRAO) and several other US universities. It primarily uses 300-m Arecibo radio 
telescope and 100-m Green Bank Telescope(GBT) to time 35 MSPs once every 
twenty days at two out of 327/400/800/1400/2350 MHz wavebands. The two fre- 
quency observations are not simultaneous, but are carried out with a gap ranging 
from 1 hour to 7 days. Both Arecibo and GBT are very sensitive telescope, thus pro- 
viding high S/N data in 15 to 45 minutes observations. Arecibo typically observes 
19 pulsars and 18 pulsars are monitored by GBT. The whole ensemble is covered 
in 18 hours per epoch. Wide-band backends, such as Astronomical Signal Processor 
(ASP) or Green Bank Astronomical Signal Processor (GASP), are used to provide 
full polar coherently dedispersed data. See Ref. f3"31 and references therein for details 
of the observing system, instrumentation used and timing analysis. 

5.3. The European Pulsar Timing Array(EPTA) 

EPTA is the most recently started PTA program. It is a collaboration between 
several European institutions, such as Max-Planck Institute for Radio Astronomy, 
Jodrell Bank Centre for Astrophysics, ASTRON and so on. It uses four major 
European radio telescopes — 70-m Lovell telescope at Jodrell Bank, 100-m Effelsberg 
telescope at Bonn, Westerbok Synthesis Radio Telescope and Nancy radio telescope. 
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Sardinia Radio Telescope is likely to be added in future. There are plans to combine 
all these telescopes coherently to form a Large European Array for Pulsars (LEAP), 
which is likely to provide a telescope equivalent to 200-m single dish. See Re fs. [5711581 
and references therein for details. 

5.4. International Pulsar Timing Array 

PTA requires frequent monitoring of pulsars at multiple radio frequencies. This adds 
up to a large observing time requirement from the participating radio telescopes. 
The pulsar community has therefore come together, in a unique initiative, to divide 
the required observing time across all participating telescopes and share the data. 



a collaboration of the three existing PTAs - PPTA, NANOgrav and EPTA, and 
encompasses almost the entire pulsar community. 

Apart from the division of observing time, such an international effort has other 
advantages and spin-offs. Firstly, this allows a fairly uniform coverage of MSPs 
in the north and south hemispheres. Secondly, the participating PTAs will have 
overlaps in their sample of MSPs, which provides consistency and quality checks on 
the data. Thirdly, the data from all these telescopes are brought together in public 
domain, which can lead to high quality data sets after proper assessment. Fourthly, 
this allows developing, comparing and contrasting alternative analysis techniques 
for data analysis through regular public data challenges. Lastly, such an effort will 
be very useful in developing the gravitational wave community. 

IPTA monitors an ensemble of about 30 MSPs with a view to detect SGWB, 
mainly from super massive black hole binaries, with a detector arm length of light 
years in the nanohertz frequency range. Thus, it is complementary to ground (space) 
based detectors with typical arm length of few (millions) km and much higher 
frequency range. 

6. Current Status of GW detection by PTAs 

The currently operating PTAs are on the threshold of detection of SGWB and are 
already putting stringent limits to the expected strain from theoretical calculations. 
A large part of PTA effort has been expended in finding the "ideal" clocks for this 
challenging and difficult experiment. With more than 20000 TOAs on all PTA pul- 
sars, timing residuals of the order of 30 ns have been obtained for PSRs J1713+0747 
and J1909-3744p2l which show almost "white" timing residuals, demonstrating that 
such high timing precision is in principle achievable. These results from NANOGrav 
are consistent with those from PPTAp^which also lists PSR J0437— 4715 as a useful 
clock. 

As discussed in Section 12. 3[ upper limits on characteristic strain, h c , can be 
estimated from these individual measurements. These results from the PTAs are 



e More information can be found at http://www.ipta4gw.org, 



This consortium is called International Pulsar 
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shown in Table [TJ It is evident that the limits have improved substantially and are 
already comparable with the expected characteristic strain due to such SGWB. 



Table 1. Summary of limits on h c from esti- 
mates on individual pulsars, assuming a=— 2/3 



Source 


h c (per year) 


Reference 


B1133+16 


< 9 X ltr 13 


Ref.QU 


B1237+25 






B1604-00 






B2045-16 






B1855+09 


< 2 x ltr 14 


Ref-ITSl 


B1937+21 






B1855+09 


< i.i x icr 14 


Ref.rgQl 


J1713+0747 


< i.i x icr 14 


Ref.1551 


J1909-3744 


< 3.9 x icr 14 


Ref.rsg 


B1855+09 


< 1.3 x icr 14 


Ref.1551 


J0437-4715 


<icr 15 


Ref.[50l 


J1713+0747 


<i(r 14 


Ref.[50l 


J1909-3744 


<icr 14 


Ref.l50l 



As discussed in Section 12.31 upper limits can also be obtained by examining 
the correlations in the timing residuals of several pairs of pulsars with varying 
angular separations. The best fit for the NANOGrav data to the Hellings & Downs 
correlation curve (Eq.EJ) is consistent with no detectable SGWB and translates to a 
2-0- upper limit of 7.2 x KT 15 , 4.1 x 10" 15 and 3.0 x 10" 15 assuming a = -2/3,-1 
and —7/6 respectively for different models of SGWB PS Likewise, PPTA data are 
also consistent with no detectable SGWB, but with a more conservative upper limit 
of 6.0 x 10" 15 P 

7. Conclusions and summary 

The current sensitivity of PTAs is tantalizingly close to the expected SGWB al- 
though no signature of such a background is evident from the data so far. The best 
"clocks" are PSR J0437-4715, J1713+0747 and J1909-3744. The most stringent 
limits in both methods employed in PTA are dominated by these pulsars. The real 
challenge for the PTA is in obtaining a larger sample of such ideal "clocks" sampling 
the correlation function in Fig. [T] uniformly. This motivates new surveys for finding 
new MSPs. New high-energy telescopes such as Fermi-LAT have already revealed a 
much larger population of MSPs and deeper and high time resolution radio surveys, 
such as HTRU survey, are underway in both the hemispheres. Another challenge is 
to rapidly characterise the new MSPs to be included in the future PTAs. Last, but 
not the least, a better understanding of magnetospheric physics of radio pulsars is 
required to model the variable torques on the neutron star in a bid to correct the 
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covariant "red noise" component in the residuals contributed by such torques. While 
the SGWB detection still seems some years away, these efforts will certainly help 
advance the understanding of physics of radio pulsars as well as their environments. 
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